AD  61 32  87 


*IW 

* 


WL  TDR-6>»-8 


WL 

TDR 

6b-8 


STEP  LOAD  MOVING  ON  THE  SURFACE  OF  A  HALF-SPACE 
OF  A  LOCKING  MATERIAL  -  SUBSEISMIC  CASE 


by 

Hans  H.  Bleich 

Paul  Weidlinger,  Consulting  Engineer 


TECHNICAL  DOCUMENTARY  REPORT  NO.  WL  TDR-6L-8 


i 


COPY  3T 

OF  _ 

HARO  COPY 

$.  5  H 

MICROFICHE 

ltd? 


Research  and  Technology  Division 
AIR  FORCE  WEAPONS  LABORATORY 
Air  Force  Systems  Command 
Kirtland  Air  Force  Base 
Nev  Mexico 


February  1965  DDC 

fnna^pnani 


DDC-IRA 


WL  TDR-6U-R 


STEP  LOAD  MOVING  ON  THE  SURFACE  OF  A  HALF-SPACE 
OF  A  LOCKING  MATERIAL  -  SUBSEISMIC  CASE 

by 

Hans  H.  Bleich 

Paul  Weidlinger,  Consulting  Engineer 


February  1965 


WL  tdr-64-8 


FOREWORD 


The  research  reported  here  was  performed  under  Project  5710,  Subtask  13. lM, 
Program  Element  7.60.06.01.5 ,  by  Paul  Weidlinger,  Consulting  Engineer,  777  Third 
Avenue,  New  York,  NY,  10017,  under  Contract  AF  29( 60l)-6055 ,  during  the  period 
July  1963  to  January  196U.  The  report  was  submitted  19  January  1965.  The  Air 
Force  project  monitor  was  lLt  Cornelius  J.  Higgins,  AFVL  (WLDC). 

This  report  has  been  reviewed  and  is  approved. 


CORNELIUS  J.  HIGGINS 
lLt  USAF 
Project  Officer 


THOMALmJ  LOWRY,  JR.^  ^  R.  A.  HOUSE 

Colonel  USAF  Colonel  USAF 

Chief,  Civil  Engineering  Chief,  Development  Division 

Branch 


ii 


WL  TDR-614-8 


ABSTRACT 


Considering  a  locking  material  prior  to  compaction  as  a  special  case  of  a 
nonlinearly  hardening  elastic  material,  conditions  at  a  discontinuity — a 
locking  front — are  analyzed  on  the  basis  of  three  dimensional  theory.  This 
study  leads  to  the  important  result  that  the  major  compressive  principal 
stress  at  a  locking  front  must  always  be  normal  to  the  front,  even  if  the 
front  is  not  plane.  Based  on  this  general  result,  the  effect  of  a  uniform 
step  pressure  traveling  with  subseismic  velocity  on  the  surface  of  a  half¬ 
space  is  obtained  for  the  case  of  a  locking  material  which  after  compaction 
has  elastic-Coulomb  behavior.  Such  a  material  acts  linearly  elastic  if  the 
state  of  stress  does  not  overcome  internal  friction,  but  slip  occurs  if  the 
stresses  reach  a  critical  state  defined  by  Coulomb  friction.  As  a  special 
case  the  solution  applies  also  for  a  material  which  is  linearly  elastic 
after  compaction.  The  stress,  velocity  and  acceleration  histories  due  to 
the  traveling  step  pressure  are  discussed  and  compared  to  those  in  the  one 
dimensional  case  of  a  suddenly  applied  uniform  surface  pressure. 
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I  INTRODUCTION 


A  preceding  report  [1],  considers  the  effect  of  a  step  pressure  pro 
gressing  with  velocity  V  on  the  surface  of  a  half-space,  Fig.  1,  when  the 
behavior  of  the  material  is  governed  by  internal  Coulomb  friction.  The  in¬ 
vestigation  [lj  studies  steady-state  solutions,  assuming  that  the  material 
acts  linearly  elastic  if  the  state  of  stress  does  not  overcome  internal  fric¬ 
tion,  but  slip  occurs  if  the  stresses  reach  a  critical  state  defined  by  inter¬ 
nal  Coulomb  friction.  Approximations  made  in  the  analysis  permit  closed  so¬ 
lutions,  but  restrict  the  result  to  velocities  V  below  a  certain  limit. 

The  present  investigation  concerns  the  equivalent  steady-state  problem, 

Fig.  2,  for  an  ideal  locking  material  which,  after  compaction  (i.e.,  after  locking) 
has  the  elastic -Coulomb  properties  assumed  in  fl|.  As  the  first  step  in  the 
analysis  a  study  of  the  locking  front  is  made.  Considering  the  locking  material 
prior  to  compaction  as  a  special  case  of  a  nonlinearly  elastic  material,  available 
three  dimensional  theory  for  nonlinearly  elastic  media  is  used  to  obtain  a  condition 
on  the  state  of  stress  at  a  front  of  discontinuity  in  such  a  medium.  Thereafter, 
the  condition  obtained  at  the  front,  and  the  approximate  procedure  developed  in 
[1]  for  materials  subject  to  Coulomb  friction  and  slip  are  jointly  applied, 
leading  to  results  in  closed  form.  The  approximation  restricts  the  results  to 
velocities  V  which  are  less  than  a  certain  limit. 

If  the  value  of  the  internal  friction  is  sufficiently  large,  slip  will  not 
occur,  and  the  present  analysis  contains  the  solution  for  the  problem  shown  in 
Fig.  2  for  a  locking  material  which  acts  linearly  elastic  after  compaction. 

The  problem  of  one  dimensional  wave  propagation  in  such  a  material  has  been 
treated  [ 2 J ,  but  the  two  dimensional  case,  Fig.  2,  has  not  been  considered  pre¬ 
viously. 
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II  THREE  DIMENSIONAL  CONSIDERATION  OF  LOCKING  FRONTS 

The  problem  of  determining  the  effect  of  a  uniform  pressure 

p(t)  acting  on  the  surface  of  a  half -space,  Fig.  3>  involves  only 

one  dimensional  strain,  so  that  it  is  sufficient  to  know  the  stress- 

strain  diagram  relating  cr  and  c  •  For  soils,  when  such  diagram 

has  the  character  of  Fig.  4a,  approximate  results  have  been  obtained 

by  replacing  the  actual  diagram  by  that  of  ideal  locking  materials, 

Fig.  4b  or  c,  [3].  The  simplification  implies  that  all  stresses  vanish 

until  the  strain  c  *  which  in  this  case  equals  the  volumetric  strain, 

XX 

has  reached  a  certain  value  defining  the  point  of  compaction. 

In  the  case  of  one  dimensional  wave  propagation,  Fig.  3>  symmetry 

requires  that  the  stress  at  the  plane  locking  front  be  normal  to  this 

front  (i.e.,  there  can  be  no  shear  stress).  The  velocity  V  of  the 

front  as  a  function  of  a  can  then  be  obtained  from  the  Rankine- 

xx 

Hugonlct  relations,  using  the  stress-strain  diagram,  Fig.  4. 

In  more  general  problems  in  two  or  three  dimensions  a  sufficiently 

small  element  of  the  locking  front  may  still  be  considered  plane,  Fig.  5, 

but  the  argument  of  symmetry  cannot  be  used  to  draw  the  conclusions  that 

the  resultant  of  the  stress  just  behind  the  front  is  normal  to  the  front 

and  that  the  shear  stress  t„  vanishes.  Nevertheless,  it  will  be  shown 

11 

below  that  the  above  conclusion  is  valid  in  the  multidimensional  case 
too,  but  the  demonstration  requires  a  much  more  involved  reasoning. 
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a)  Non-linear  Elastic  Materials 


Postponing  the  complications  due  to  the  effect  of  slip  in  a  Coulomb 
type  material,  this  section  is  restricted  to  non-linear  elastic  materials 
which  harden  with  reduction  in  volume.  Conclusions  drawn  can  then  be 
applied  to  ideal  locking  materials,  because  they  are  special  cases  of  such 
a  non-linearity. 

To  study  locking  fronts,  a  mathematical  description  of  three  dimensional 

behavior  is  required  for  a  homogeneous,  non-linear  material  which  acts  in 

one  dimensional  strain  according  to  Fig.  4a.  Using  the  concept  of  a  ''strain- 

energy  function"  W,  [4,  p.  95J,  the  stresses  o.  in  an  elastic  material  may 

1  J 

be  obtained  as  derivatives  of  W  with  respect  to  the  strains  , 


aw  aw 

°ii  =  a«ijL  '  aij  =  sat.j 


(2-1) 


where  the  subscripts  i  and  j  become,  in  turn  the  Cartesian  coordinates 
* 

x,  y  and  z.  For  isotropic  materials  the  strain  energy  function  must  be 
invariant  to  transformation  of  the  coordinates,  which  requires  W  to  be  a 
function  of  the  invariants  of  the  strain  tensor.  For  linear  elastic 
materials  this  leads  [4.  p.  102]  to  the  function 


Note  that  [4]  uses  different  symbols  and  a  different  definition  of  the 

shear  strain.  To  obtain  the  equations  in  this  paper  the  following 

substitutions  were  made:  The  stresses  in  [4],  X  ,  X  ,  etc.,  become 

x  y 

0  ,0  ,  etc.,  respectively;  longitudinal  strains  e  ,  e  ,  e 

xx  xy  xx'yy'zz 

become  «  ,  «  ,  e  ,  while  the  shear  strains  e  ,  e  ,  e  become 

xx  '  yy  ’  zz  xy  7  yz  zx 

2t  ,  2c  ,  and  2c  .  As  a  result  of  the  change  in  definition  of 
xy  *  yz  zx 

the  shear  strains,  the  quantities  e, ,  are  components  of  a  tensor. 

■4 1, 


3 


(2-2) 


w  =  ( i x)2 


2n  I2 


where  K  and  p  are  constants,  while 


e  +  e  +6 
xx  yy  zz 


c  e  +  e 
yy  zz  zz 


e  +  e 

XX  XX 


2 

e 

xy 


(2-3) 


1^  and  Ig  are  invariants  of  the  strain  tensor.  It  is  noted  that  1^ 
describes  the  change  in  volume  of  an  element. 

To  describe  the  most  general  case  of  an  isotropic  non-linear 
material,  the  constants  X  and  n  in  Eq.  (2-2)  should  be  considered  as 
functions  of  all  invariants  of  the  strain,  of  which  there  are  three. 
In  addition  to  1^  and  Ig,  given  by  Eqs.  (2-3)>  there  is  a  third  one 
[5,  p.  180,  Eq.  (25.11)] 


e  e  e  +  2c  e 
xx  yy  zz  xy  yz 

2  2 
Gyy  Sxz  “  ezz  exy 


zx 


e 


XX 


(2- 3a) 


However,  in  order  to  describe  a  material  which  becomes  harder  with 
consolidation,  i.e.,  with  volume  change,  \  and  M-  should  be  just  func¬ 
tions  of  1^,  which  defines  the  volumetric  strain.  It  will  be  shown 
below  that  for  the  present  purpose,  one  may  assume  that  K  and  n  are 
proportional, 


k 


x  =  x0  HIX) 


^  =  %  ^i) 


(2-4) 


where  X  and  are  positive  values  having  the  dimensions  of  Lame's 
constants  X  and  \i ,  and  the  function  F(l^)  is  always  positive.  In 
order  to  represent  a  hardening  material,  the  derivative  of  F  must 
satisfy  the  inequality 


*  0 


(2-5) 


when  <  0.  Substitution  of  Eqs.  (2-4)  into  (2-2)  gives 


W 


X  +2u  5 

-V^  p(Ii>  h  -  2(1o  h 


(2-6) 


Application  of  Eqs.  (2-1)  furnishes,  after  rearrangement,  the  stress- 
strain  relations 

“ii  “  h  +  2|io  eii>  +  I  +  ^  h2  -  K  I2] 

(2-7) 

and  for  i  ^  j 


where 


°ij  "  2F(I1)  ^o  ei(j 


dF(lJ 


F'^I1^  “  dl. 


(2-0) 


(2-9) 
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To  demonstrate  that  the  assumption  made  in  writing  Eqs.  (2-4) 

is  permissible,  it  will  be  shown  that  it  leads  to  a  diagram  of  the 

type  of  Fig.  4  in  the  case  of  uni-axial  strain  in  the  direction  of 

the  x-axis  when  all  values  e .  .  vanish  except  e  .  For  this  case 

ij  xx 

Eqs.  (2-3)  give  I1  =  and  1 2  =  0,  such  that  F(l1)  =  F(sjcx)> 

F'(I, )  =  F'(e  ).  The  stress-strain  relation  is  obtained  by  substi- 

_L  XX 

tuting  these  expressions  into  Eq.  (2-7)  for  i  =  x: 


xx 


=  (X  +  2p,  ) [e  F(e  )  +  J  e2  F'(e  )] 
'  o  xx  '  xx'  2  xx  '  xx' 


(2-10) 


By  definition,  the  function  F  is  positive  and  by  virtue  of  Eq.  (2-5) 
F  increases  for  negative  (compressive)  strains  with  j  e  l  .  The  term 
e  F(e  )  in  Eq.  (2-10)  has  therefore  the  character  of  Fig.  4a,  in- 

XX  XX 

eluding  the  negative  sign  in  the  range  e  <  0.  The  other  term, 

1  2 

-  e  F'(e  )  is  also  negative.  To  assure  that  its  absolute  value 
c  XX  XX 

also  increases  similarly,  we  prescribe 

F"(e  )  >  0  (2-11) 

XX 

Both  terms  having  the  character  of  Fig.  4a,  the  entire  expression 
will  have  the  same  character.  The  strain-energy  function  (2-6)  fur¬ 
nishes  therefore  a  three  dimensional  expression  for  materials  having 
a  one  dimensional  stress-strain  diagram  of  the  type  of  Fig.  4,  pro¬ 
vided  the  function  F  satisfies  Eqs.  (2-5)  and  (ll).  [An  example  of 
a  function  satisfying  the  requirements  is  F  =  exp(-Ne  ),  where  N 
is  a  positive  number.  It  is  easily  verified  that  Eq.  (2-10)  gives 


6 


a  a-t  diagram  of  the  type  of  Fig.  4a. , 


The  stress-strain  relation  (2-1)  and  the  expression  (2-6)  for 
W  can  now  be  introduced  into  the  Rankine-Hugoniot  relations.  Con¬ 
sider  two  dimensional  wave  propagation  for  the  case  of  plane  strain. 
Let  the  x  and  y  axes,  respectively,  be  normal  and  tangential  to  the 
front,  Fig.  6,  while  the  z-axis  is  normal  to  the  plane  of  the  paper. 


Fig.  6  indicates  a  rectangular,  infinitesimally  small  element, 

in  its  original  shape  ABCD,  and  in  its  distorted  shape  A'B'CD,  after 

the  wave  front  has  passed  over  it.  Let  U  be  the  propagation  velocity 

of  the  front.  In  the  time  dt  the  wave  front  has  traveled  the  distance 

U  dt  =  AD,  while  point  A  has  traveled  to  A’,  the  x  and  y  components 

of  its  motion  being  u  dt  and  u  dt,  respectively,  where  u  ,  u  are 

^  y  x  y 

the  components  of  the  particle  velocities.  From  Fig.  6  one  may  read 
directly  the  longitudinal  strain 


xx 


u 


U 


(2-12) 


The  minus  sign  is  due  to  the  fact  that  the  positive  directions  of 
a  and  u  are  opposite  each  other.  Noting  further  that  the  con- 
ventional  shear  strain  y  equals  2e^.  ,  one  finds  similarly 

u 

7  =2*  =  - 
'xy  Cfxy  U 


.(2-13) 


In  the  limit  dt  -*  0  the  distance  AB  equals  A'B'  such  that  on  the  shock 

front  c  =0.  Further,  because  of  plane  strain,  «  =  e  =  e  =0 

yy  *xz  *yz  *zz 
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such  that  only  e  .  and  e  need  be  considered  further,  Eqs.  (2-3) 
xx  xy 

become  therefore 


I 


2 


(2-14) 


and  from  Eq.  (2-7): 


a  =  (X  +  2n  )  e  P(  e  )  +  i  [  (X  +  2p.  )  e  2  +  4|x  e  2  ]  F '  ( e  ) 
xx  o  o  xx  v  xx'  2  v  o  o'  xx  o  xy  xx 

(2-15) 

a  =■■  a  =  X  e  F(e  )  +  4  [(X  +  2u  )  e  2  +  4u  e  2  ]  F'(e  ) 

yy  zz  o  xx  v  xx;  2  v  o  ^o'  xx  o  xyJ  v  xx' 

(2-16) 

The  shear  stresses  are  obtained  from  Eq.  (2-8): 


a  =  2|i  {  F(e  ) 
xy  o  xy  '  xx' 


(2-17) 


a  =  a  =0 
xz  yz 


(2-18) 


In  the  time  dt  in  which  the  front  engulfs  the  element  ABCD, 

Fig.  6 ,  the  x  and  y-conrponents  of  the  velocities  of  all  the  particles 

in  the  element  are  changed  by  u  and  u  ,  respectively.  The  mass  of 

x  y 

the  element  being  p  U  dt  dy  dz,  the  change  in  momentum  in  the  x- 
direction  can  be  expressed  by  the  stress  a  , 

Aa 


p  U  u 

^  X 


XX 


(2-19) 
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where  p  is  the  density.  A  similar  relation  applies  for  the  y- 
direction 


p  U  u  =  -  a 
1 


xy 


(2-20) 


The  minus  signs  appear  again  because  the  directions  of  the  stresses 
shown  in  Fig.  6  are  opposite  to  the  positive  directions  of  the  velo¬ 
cities.  Substituting  the  values  of  u  and  u  into  Eqs.  (2-12,  13), 

^  y 

one  finds 


o  =  p  u  e 
xx  xx 


u2 


(2-21) 


a  =  2p  IT  e 
xy  xy 


(2-22) 


Eqs.  (2-21)  and  (2-15)  give 


F1  (e  ) 

U2  =  (X  +  2ul  )  F(e  )  +  i  [(X  +  2[i  )  c2  +  ka  e  2  ]  - — 

'  «  '  vv'  2  '  o  o'  xx  o  xy  e 


O  O'  XX' 


XX 

(2-23) 


The  division  by  e  ,  required  to  obtain  this  result,  is  always  per- 

mitted,  because  a  consolidation  front  is  being  investigated,  such 

that  e  ^  0.  Substituting  Eq.  (2-22)  into  (2-17)  gives,  after 

division  by  2s  ,  the  alternative  relation 
xy' 


P 


F(e  ) 
'  xx 


(2-24) 


This  equation  applies  only  if  e 


*y 


0,  the  division  not  being  possi- 
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ble  if  e  =0. 
xy 

If  the  shear  strain  e  does  not  vanish,  Eqs.  ('2-23,  24) 

xy 

should  both  be  valid;  however,  it  can  be  shown  easily  that  this  is 
not  the  case.  Noting  Eq,  (2-5)  and  <  0,  XQ  >  0,  >  0,  both 

terms  on  the  right  hand  side  of  Eq.  (2-23)  are  positive,  such  that 

P  U2>  (XQ  +  2v.Q)  F(exx)  (2-25) 

This  is  incompatible  with  Eq.  (2-24).  The  contradiction  disappears 

only  if  e  =0.  In  this  case  Eqs.  (2-22  and  20 )  give 
xy 

c  =  0  ,  u  =  0  (2-26) 

xy  '  y 


It  is  therefore  concluded  that  on  a  compaction  front  the  stress  re¬ 
sultant  and  the  particle  velocity  are  normal  to  the  front. 

Introduction  of  e  =0  into  Eqs.  (2-15,  16)  gives  the  follow- 
xy 

ing  expressions  for  the  stresses  at  the  front: 


xx 


\  +2[i 

=  (\  +  2p.  )  e  F(e  )  +  -°  0 


xx 


XX' 


e  F'(e  ) 
xx  xx' 


(2-27 ) 


yy 


=  cr. 


zz 


= 


XX 


F(e  ) 
v  xx' 


e  F'(e  ) 
xx  '  xx' 


(2-28) 


Xq,  and  F  being  by  definition  positive,  while  <  0  (compac¬ 
tion)  smd  F‘  <  0  [see  Eq.  (2-5)3,  all  terms  on  the  right-hand  side 
of  Eqs.  (2-27,  28)  are  necessarily  negative.  Using  these  facts, 


10 


I 


the  relative  magnitudes  of  a  and  a  =  a  „  may  be  compared.  The 

xx  yy  22  ^ 

expressions  for  these  stresses  differ  only  in  the  term  containing 

F(e  ),  and  one  finds  the  inequality 
XX  * 


b)  Locking  Materials ,  Llnearily  Elastic  after  Compaction 


The  results  of  Section  a)  will  now  be  applied  to  an  ideal 
locking  material  having  a  uni-axial  stress-strain  diagram  of  the 
type  Pig.  4c,  i.e.  a  material  which  becomes  linearily  elastic  after 
locking.  Using  a  limiting  process,  the  diagram  4c  is  being  considered 
as  a  special  case  of  the  diagram  in  Fig.  7>  in  the  limit  The 

XX 

material  corresponding  to  Fig.  7  Is  defined,  as  follows: 

For  compressive  strains  between  0  and  -c,  the  stresses  are 

obtained  from  the  noi.  -linear  relations  (2-27,  28),  reaching  values 

oL  and  o ^  for  e  =  -e.  For  strains  beyond  this  limit  the 

xx  yy  zz  xx 

material  is  deemed  to  be  linear,  having  material  constants  E  and  v. 

The  stress  increments  Acr^  after  locking, 


Aoii  =  aii  ~  aii 


(2-30) 


can  be  obtained  from  the  strain  increments  Ae^  beyond  the  locking 
strain  by  using  the  conventional  elastic  relations.  The  x-component 
of  the  strain  change  is 


Ae  =  e  +  e 

XX  XX 


(2-31) 


All  other  components  of  the  strain  vanish  at  the  front,  Eqs. 
(2-13a,  13b,  26a),  such  that  all  increments  Ae.^  also  vanish,  except 
Ae  .  In  this  situation  (transverse  strains  being  inhibited),  the 

XX 

f 

ratio  of  the  transverse  and  longitudinal  stresses  according  to  linear 

v 

theory  of  elasticity  is  ,  such  that 
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t 


ioyy  -  iaZZ  *  Cv  ioxx 


(2-38) 


Consider  now  Eq.  (2- 30)  in  the  limit 


lim  a  -40 
xx 


(2-33) 


Because  the  inequality  (2-29)  requires  the  absolute  values  of  a. 


yy 


and  a  to  be  less  than  a  ,  it  is  necessary  that 
zz  xx 


lim  a L  =  lim  0 
yy  zz 


(2-33b) 


Going  to  the  limit  0,  Eq.  (2-30 )  becomes  therefore 


°ii  -  4o,ii 


(2-3*1-) 


and  substitution  into  Eq.  (2-32)  leads  to  a  relation  for  the  stresses 
just  behind  the  consolidation  front 


a  *>  a  =  a  (2-35) 

yy  zz  1-v  xx  v  ' 


while  according  to  Section  a)  all  shear  stresses  vanish: 


a  =  a  =  a  =0  (2-35a) 

xy  yz  zx  '  ' 


Tliese  relations  apply  immediately  behind  the  compaction  front,  and 
will  be  used  in  Section  III  as  boundary  conditions  when  determining 
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the  stresses  after  compaction  at  points  anywhere  behind  the  front. 
The  value  of  the  increment  of  the  strain  Ae  in  terms  of 

XX 

A a  obtained  from  the  conventional  elastic  relations  is 

XA 


2 

1-V-2V* 

E(l-v) 


A  a. 


xx 


(2-36) 


or  from  Eqs.  (2-31)  and  (2-34) 


l-v-2v  ‘ 

exx  e  +  E(l-v)  axx 


(2-37) 


This  relation  applies  for  the  stress-strain  diagram  Fig.  4c.  To 
obtain  the  relation  corresponding  to  Fig.  4b.  let  E-»a>,  with  the  re¬ 
sult 


S*  -  '  e  (2-38) 

Tills  equation  defines  the  strain  immediately  behind  the  compaction 
front.  Substituting  this  value  into  Eq.  (2-21)  gives  the  velocity 
U  of  the  front 


U  - 


(2-39) 


Due  to  a  <  0  (compression),  the  expression  under  the  root  is  posl- 
tlve,  as  it  should  be.  (Note  that  only  the  positive  value  of  the 
square  root  has  physical  meaning.  ) 


Ik 


c)  Locking  Materials,  Linearly  Elastic  but  Subject  to  Slip  after 
Compaction 

Consider  finally  locking  fronts  in  two  dimensional  problems 
of  plane  strain  for  an  ideal  locking  material  which  behaves  linearly 
elastic  provided  internal  friction  of  the  Coulomb  type  is  not  over¬ 
come,  but  which  will  slip  when  the  stresses  satisfy  Coulomb’s  rule. 

Ref.  [1,  p.  6,  73  contains  for  this  case  a  formulation  of  the 
slip  condition  which  is  convenient  for  the  present  study  of  the  state 
of  stress  at  the  compaction  front.  Let  o^,  Og  be  the  major  and  minor 
principal  stresses  in  the  plane  of  action,  both  being  restricted  to 
compression,  such  that 

Of  <  °  ,  a2<  0  ,  -ff1  *  -ag  (2-40) 

Defining  a  quantity  a  by  the  relation 

a  =  -  (o^  +  s)  (2-41) 

the  possible  states  of  stress  in  the  material  are  further  restricted 
to 

-Og  £  k  c  (2-42) 

where  k  is  a  positive  number,  k  <  1,  which  is  related  to  the  angle 
of  interior  friction,  while  s  >  0  defines  the  value  of  cohesion. 

No  slip  can  occur  if 

-o2  >  k  a  (2-43) 
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while  slip  occurs  if 


-0g  =  k  a  (2-44) 

Values  of  -cr^  <  k  a  can,  by  definition,  not  occur  in  this  material. 
The  above  conditions  involve  only  the  principal  stresses 
and  a^,  in  the  plane  of  the  tiro  dimensional  problem.  In  order 
that  this  formulation  remains  valid  the  third  stress  a ^  lie  te“ 

tween  the  other  two,  -c^  ^  -cr^  ^  ~a2‘ 

Concerning  the  strains,  it  was  assumed  in  [l]  that  the  plas¬ 
tic  (or  slip)  deformations  occur  without  volume  change,  while  the 
elastic  deformations  are  negligible.  In  the  present  case,  a  more 
refined  definition  must  be  derived,  allowing  for  the  elastic  strains. 

Retaining  the  result  obtained,  in  Section  a),  stating  that  the 
stress  resultant  at  the  compaction  front  is  normal  to  this  front, 
the  state  of  stress  and  the  associated  strains  can  be  obtained  re¬ 
latively  simply.  Consider  a  rectangular  element  ABCD,  Fig.  8,  as 
the  compaction  front  passes  over  it,  and  let  x  be  a  coordinate  at 

right  angles  to  the  front.  On  the  front,  the  tangential  strains 

T 

must  vanish  such  that  the  only  non- vanishing  strain  is  e  ,  where 

AA 

the  superscript  T  indicates  total  strain  from  all  causes.  For  an 
elastic-plastic  material  this  total,  strain  is  the  sum  of  an  elastic 
contribution  (in  the  present  case  non-linear)  and  of  the  slip  strain, 
and  we  have  the  relations 
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e  +  e  ® 

XX  XX 


q  m 

e  +  e  =e  ==0 
zz  zz  zz 


(2-45) 


g 

where  are  the  elastic  strains,  the  slip  strains,  and  y,z  are 
coordinates  in  the  plane  of  the  front.  Adding  the  three  equations 
gives 


=  (e  +  e  +s  )  +  (e°+e^+e^)  (2-46) 

xx  '  xx  yy  zz'  v  xx  yy  zz'  v  ' 


The  sum  of  the  three  strains  e. .  is  the  volumetric  strain  e  ,  and 

ii  v7 

g 

similarly  the  sum  of  the  three  strains  ei;L  is  the  volume  change  dur¬ 
ing  slip.  In  conformity  with  [l]  the  latter  volume  change  vanishes, 
such  that 


(2-47) 


T 

The  total  strain  e  .  is  therefore  equal  to  the  elastic  volume  change. 

XX 

For  an  ideal  locking  material,  corresponding  to  Fig.  4b,  the  volume 
change  does  not  depend  on  the  stress  level,  =  -e,  or 


(2-48) 


where  e  is  the  locking  strain.  Eq.  (2-48)  is  identical  with  Eq. 
(2-38)  and  the  velocity  of  the  compaction  front  is  again  given  by 
Eq.  (2-39). 
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The  state  of  stress  in  case  of  slip  is  of  course  different 

from  the  elastic  case,  Eq.  (2-35).  Using  Eqs.  (2-4l  and  44),  the 

major  principal  stress  is  a,  ==  a  ,  while  a0  =  a  ,  and  one  finds 

jl  xx  ^  yy 

the  tangential  stress 


a  -  k(cr  +  s) 
yy  v  xx  ' 

If  cohesion  vanishes,  s  =  o,  and 


(2-49) 


k  a 


xx 


(2-50) 


This  value  of  the  stress  a  applies  immediately  behind  the  com- 

yy 

paction  front,  and  will  be  used  later  as  boundary  condition  when 

determining  the  stress  field  behind  the  front. 

Eqs.  (2-49)  and  (2-50)  apply  only  if  slip  actually  occurs, 

which  depends  on  the  relative  values  of  k  and  v.  To  decide  if  slip 

occurs,  consider  Eq.  (2-35)  for  a  without  slip, 

*y  «y 


V 

a  =  s —  a 
yy  1-v  xx 

If  in  the  case  without  cohesion 


(2-51) 


*  k  (2~52) 

Eq.  (2-43)  is  not  satisfied,  such  that  slip  must  occur  and  Eq.  (2-50) 
is  to  be  used.  If  on  the  other  hand 
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(2-53) 


v 

1-v 


>  k 


Eq.  (2-43)  is  satisfied  such  that  the  behavior  a.  the  front  is  elas¬ 
tic  and  the  elastic  relation  (2-35)  applies. 
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III  EFFECT  OF  A  PROGRESSING  STEP  LOAD  ON  THE  SURFACE  OF  A  HALF-SPACE 
The  principal  subject  of  this  study  is  the  case  of  locking  materials  which 

i 

?  may  slip  after  compaction.  However,  for  purposes  of  presentation  it  is  more 

I 

;  convenient  to  consider  first  the  less  interesting  case  of  a  material  which, 

1 

‘  after  locking,  acts  elastically,  without  slip.  The  latter  material  is  con- 

j  sidered  in  this  Section  as  Case  A,  while  the  material  with  slip  is  considered 

I 

as  Case  B. 

Case  A.  After  Locking, the  Material  is  Elastic. 

Figure  9  indicates  the  progressing  step  load  and  a  system  of  polar  coor¬ 
dinates  9,  r  moving  with  the  load.  Dimensional  considerations,  similar  to  those 
used  in  [1],  require  that  in  the  present  problem  the  stresses  be  proportional 
to  the  load  p,  but  otherwise  solely  functions  of  the  variable  9, 

o,  t  ...  =  p  f(cp)  . 

The  fact  that  the  stresses  are  only  a  function  of  9  implies  that  the  locking 
front  is  a  plane,  inclined  at  an  as  yet  unknown  angle  9=0. 

The  state  of  stress  in  the  wedge  of  angle  0  behind  the  front  is  governed 

by  the  appropriate  equations  of  elasticity.  In  similarity  to  [1],  the  present 

analysis  will  be  restricted  to  velocities  V  which  are  a  fraction  of  the  velocities 

c  and  c  of  wave  propagation  in  the  elastic  material  after  locking.  The  differ- 
P  s 

ential  equations  for  the  stresses  in  polar  coordinates  for  this  case  are  given 
in  [1,  Sq.  (A-4,5),  p.  42  and  43] 
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(3-1) 


dt 

dxp 


+  o 


r 


0 


da 

J  + 

ckp 


2t 


0 


+  a  )  *  0 
<P 


where  a  ,  a  are  the  normal  stresses 
r  ’  cp 

the  shear  stress.  On  the  surface  cp  * 


in  the  r-  and  9-  directions,  while  t  is 
0,  the  boundary  conditions  are 


a  (0)  -  -  P  7(0)  -  0  (3-2) 

while  behind  the  compaction  front,  at  cp  *  8,  the  conditions  (2-35)  obtained 
in  Section  lib  apply.  Noting  that  the  coordinates  x,  y,  z  in  these  conditions 
are  normal,  and  parallel  to  the  front,  respectively,  the  stress  symbols  are  to 
be  interpreted  as  follows: 


a  =a  ,  a  =  a  ,  a  s  t 
XX  cp  yy  r  xy 


such  that  one  has  two  additional  conditions: 


"  (3) 


v 


r'r  *  1-v  9 

t(8 )  -  0  . 


o  (3) 


(3-3a) 

(3~3b) 


The  solution  of  the  differential  equations  (3-1)  satisfying  the  two 
conditions  (3-2)  may  be  taken  directly  from  Ref.  [1]  (Eq.  A-7), 

c^  =  -  p  -  D(1  +  cos  29)  -  0(29  +  sin  29) 

a  «  -  p  -  D(1  -  cos  29)  -  0(29  -  sin  29)  (3_iO 

T  «  D  sin  29  +  C(l  -  cos  29) 
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where  C  and  D  are  constants  to  be  determined  from  Eqs.  (3-3)*  It  can  be 
verified  by  substitution  that  the  second  of  these  equations  is  satisfied  if 

C  »  -  N  cos  3  D  «  N  sin  3  (3~5) 

where  N  is  a  new  constant  to  be  determined  from  the  remaining  boundary 
condition  (3-3a).  Substitution  of  Eqs.  (3-5)  and  (3-3a)  into  Eqs.  (3-*0> 
gives  two  simultaneous  linear  equations  on  N  and  0^(3)*  Determining  0^(3) 
one  finds 


-  (1-v)  P _ 

v  +  (1-2 v)  (3  cot  3 


(3-6) 


This  equation  contains  the  still  unknown  angle  3*  To  determine  the  latter 

the  relation  (2-39)  for  the  velocity  U  of  the  shock  and  compaction  front 

is  to  be  used  (note  a  =  o  ) 
v  xx  cp ' 


P 


(3-7) 


The  velocity  U  of  the  front  must  be  the  component  at  a  right  angle  to  this 
front  of  the  velocity  V  of  the  load,  requiring 


U  -  V  sin  3 


(3-8) 


Eqs.  (3-6,  7 )  and  8)  give 

V2  -  -2-  - -p - 1  -  (3-9) 

pe  v  sin  3  +  (l-2v)  sin  3  cos  3 

This  equation  defines  3  in  terms  of  V.  It  is  simpler,  however,  to  use  this 
equation  to  compute  V  for  given  values  of  3*  If  the  relationship  between  V 
and  3  is  plotted  for  meaningful  values  of  0  <  v  <  £  ,  a  curve  of  the  type  shown 
in  Figure  10  for  v  »  l/3  is  obtained.  The  velocity  V  has  a  minimum,  Vcr  ,  and 


22 


there  are  two  values  of  the  angle  3  for  each  value  V  >  V  .  The  critical 

cr 

minimum  value  of  V  occurs  for  a  critical  value  of  the  wedge  angle  3cr  <  tt/2. 

p 

(There  are  also  values  3  >  2.03  in  Fig.  10  where  V  <  0.  This  range  of 
imaginary  velocities  V  is  physically  meaningless.) 

3  as  function  of  V  is  double -valued,  but  can  be  made  single-valued  through 
elimination  of  one  of  the  values  3  by  qualitative  reasoning.  One  is  interested 
in  a  steady  state  which  is  the  long-time  solution  to  a  load  p  applied  to  an 
expanding  region,  such  that  3  should  be  less  than  tt/2,  Figure  11.  If  3  >  tt/2, 

Figure  12,  one  would  find  the  following  unreasonable  situation.  Just  behind  the 
shock  front,  the  particle  velocity  is  necessarily  at  right  angles  to  the  front, 
and  in  case  of  Figure  12,  the  velocity  would  be  upward.  As  the  particle  velocity 
further  behind  the  front  changes  only  slightly  from  the  value  at  the  front,  the 
downward  loads  p  would  produce  upward  velocities,  which  is  not  possible.  Among 
the  steady-state  solutions  contained  in  Eq.  (3-9)>  those  for  3  >  tt/2,  while 
mathematically  correct,  are  physically  not  appropriate.  By  reasoning  further,  that, 
in  the  range  ®  >  V  >  where  solutions  exist  at  all,  the  solutions  ought  to  be 
continuous,  we  conclude  that  the  solutions  for  angles  in  the  range  3cr  <  3  <  tt/2 
should  also  be  eliminated,  and  only  the  range 

0  <  3  <  3cr  (3-10) 

furnishes  meaningful  solutions.  In  the  example,  Figure  10,  the  applicable 

values  of  V  are  shown  as  a  solid  line.  The  value  of  p  can  be  found  by  aifferen- 

cr 

tiation  of  Eq.  (3-9 ) >  which  leads  to  the  condition 


tan  (20cr) 


-  2(l-2v) 


(3-H) 


The  fact  that  no  steady- state  solution  was  found  if  the  velocity  V  is 
less  than  a  critical  value,  Vcr  ,  requires  discussion.  If  one  considers  the 
problem  of  a  half-space  subjected  to  an  expanding  load  p,  the  solution  may  have 
the  character  shown  in  Figure  11,  where  the  locking  front  and  the  leading  edge 
of  the  pressure  p  meet  at  point  A.  However,  it  is  also  possible  that  the  locking 
front  is  a  "detached"  front,  Figure  13,  moving  ahead  of  the  load  p.  In  this  case 
the  "steady-state"  solution,  which  is  an  approximate  solution  near  the  front  of 
the  applied  load,  point  A,  is  simply  the  Huth-Cole  solution  [6],  because  in  the 
vicinity  of  A  the  material  is  elastic  and  the  compaction  front  has  moved  far 
away. 

Substituting  B  from  Eq.  (3-11)  into  Eq.  (3”9)>  it  will  be  noted  that 

vJL 

is  proportional  to  p/p*  such  that  the  critical  velocity  V  is  proportional  to 

OX* 

/p  ,  and  V  increases  with  the  applied  pressure, 
cr 


Equation  (3-6)  permits  the  determination  of  the  stress  at  the  locking 
front,  the  value  of  which  may  be  shown  to  be  the  maximum  of  or^  anywhere.  Figure  14 
gives  a  plot  of  -  cr ^p  as  function  of  V,  indicating  that  varies  only  slightly, 
from  unity  for  V  ■  »  t  up  to  a  -1.37  for  V  *■  V 


<P 


cr 


Case  B.  After  Locking,  the  Material  is  Elastic,  but  Subject  to  Slip. 

For  this  material,  the  conditions  at  the  compaction  front  have  been  con¬ 
sidered  in  Section  lie.  The  relations  between  the  stresses  at  the  front  depend 
on  the  relative  values  of  Poisson's  ratio  v  and  of  the  material  parameter  k. 


<  k 
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If,  see  Eq.  (2-51), 


_v_ 

1-v 


(3- 12) 


slip  will  occur  in  an  infinitely  thin  region  immediately  at  the  locking  front, 
and  the  relations  between  the  normal  and  tangential  stresses  immediately  be¬ 


hind  the  front  are  given  by  Eq.  (2-50).  If  on  the  other  hand, 

ir;  >  k  (3-13) 

the  material  at  the  front  will  remain  elastic,  and  Eq.  (2-35)  applies. 

Consider  first  materials  where  Eq.  (3-12)  holds,  such  that  Eq.  (2-50) 
applies.  Using  again  the  system  cf  polar  coordinates  r,  cp  shown  in  Fig.  9> 
the  state  of  stress  in  the  wedge-shaped  region  0  <  cp  <  0  will  be  subject  to 
the  following  boundary  conditions: 


ar(e) 

*  k  a  (3) 

Y 

(3-l4a) 

t(0) 

-  0 

(3-l4b) 

o  (0) 
<P 

»  -  p 

(3-l4c) 

t(0) 

»  0 

(3-l4d) 

In  the  interior  of  the  wedge,  the  behavior  of  the  material  was  des¬ 
cribed  in  Section  lie  by  a  set  of  equations  and  inequalities,  Eqs.  (2-40  to  44). 
To  simplify  the  present  treatment  materials  without  cohesion  will  be  considered, 
such  that  the  quantity  s  in  Eq.  (2-4l)  vanishes,  s  ■  0. 


Equations  (3-l4)  give  conditions  for  the  stresses  at  the  boundaries  of 


the  wedge  formed  by  the  surface  and  by  the  compaction  front,  but  the  determi¬ 


nation  of  the  stresses  within  the  wedge  is  complicated  by  the  fact  that  one  does 
not  know  in  this  stage  in  what  portion  of  the  wedge  slip  occurs,  and  in  what 
portions  the  material  will  act  elastically.  It  will  be  demonstrated  hereafter 
that  elastic  behavior  occurs  everywhere  in  the  interior  of  the  wedge. 
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To  show  this,  the  elastic  solution  will  be  obtained  from  the  differential 
Eqs.  (3-1 )  used  in  Section  A,  in  conjunction  with  the  boundary  conditions 
(3-1*0  appropriate  for  the  case  of  slip  at  the  compaction  front. 


Noting  that  the  present  boundary  conditions  (3-l4c,  d)  are  identical 
with  Eqs.  (3-2)  in  Section  A,  the  solution  (3-A)  of  the  differential  Eqs.  (3-1) 
applies  again, 

■  -  p  -  D(1  +  cos  2tp)  -  C(2tp  +  sin  2cp) 
a  -  -  p  -  D(1  -  cos  2cp)  -  C(2cp  -  sin  2tp)  (3-15) 

T  «  D  sin  &p  +  C(1  -  cos  2<p) 

where  the  arbitrary  constants  C  and  D  will  be  determined  from  the  boundary 
conditions  (3-l4a,  b).  The  second  of  these  conditions  is  satisfied  if 

C  ■  -  N  cos  p  D  «  N  sin  P  (3-16) 


where  N  is  a  new  constant  to  be  determined  from  Eq.  (3-l*+a),  Substitution  of 
Eqs.  (3-lta)  and  (3-16)  into  Eq.  (3-15)  for  q>  »  P  gives  two  simultaneous 


equations  on  N  and 


(b). 


One  finds 


(3-17) 


2P  C08  P  + 


sin  P 


(3-18) 


To  confirm  the  claim  that  the  stresses  everywhere  in  the  wedge  are  such 
that  no  slip  occurs  and  that  the  solution  (3-15)  is  appropriate,  it  is  con¬ 
venient  to  express  the  stresses  and  r  in  terras  of  the  principal  stresses 


26 


and  02  and  of  the  angl*  7  ■  7(cp)  which  the  direction  of  0^  makes  with  the 

,  radial,  Figure  15.  Proceeding  similar  to  the  approach  in  [1,  p.  43,  Eqs. 

1  (A-8  to  10)],  the  principal  stresses  0-  and  cf^  are  expressed  by  two  functions 

1  d 

of  cp,  namely  0  *  a(cp )  ann  k  *  k(cp), 


(3-19) 


o2  *  -  ko 


Using  the  conventional  ( xpressions  for  the  stresses  in  terms  of  the  principal 
stresses,  one  finds 

2  2  -  2 
cos  7+0^  sin  7  «  -  (l-k)o  cos  7  -  ko 


2  2  -  2  - 

0-0..  sin  7  +  0  cos  7  *  -  (l-k)o  sin  7  -  ko 

Cp  x  d 

r  ■  (0^  -  0g)  sin  7  cos  7  *  -  (l-ic)a  sin  7  cos  7 


(3-20) 


Elimination  of  0  and  1<  from  these  equations  yields  a  relation  between  , 

0  ,  t  and  7 

<P  ' 


0  -  0 

r  cp 


tan  27-0 


(3-21) 


Alternatively,  0  may  be  eliminated  from  the  first  two  Eqs.  (3-20)  with  the 


result 


±-t. 


,  0-0 

1  r  cp 


(3-22) 


1+k  cos  27  0  +  0 

r  cp 


The  stresses  ,  0^  and  7  for  any  value  of  cp  are  given  by  Eqs,  (3-15) 
I  if  the  values  of  C,  D,  and  N  from  Eqs.  (3-16,  18)  substituted.  To  show 
that  no  slip  occurs  in  the  interior  of  the  wedge,  it  will  be  &hown  that  the 
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ratio  of  the  principal  etresses  obtained  from  Eqs.  (3-19) 

o2/a1  -  k  (3-23) 

remains  within  required  limits.  If  0^  is  the  major  stress,  -  slip 

will  not  occur  if  k  >  k.  However,  the  formulation  (3-20)  does  not  assure 
that  is  the  major  principal  stress.  If  0g  is  the  major  stress,  slip  will 
not  occur  if  ic  <  l/k,  such  that  the  general  requirement  for  elastic  behavior 
without  slip  becomes 

l/k  >  k  >  k  (3-24) 

In  similarity  to  [1,  Appendix  A],  this  inequality  will  be  confirmed 
step  by  step  in  various  locations.  As  a  first  step,  the  behavior  of  the 
function  k(cp)  near  9  «  3  is  considered. 


Due  to  the  boundary  conditions  (3-l4a,  b),  and  a r  in  the  location 
9*3  are  the  major  and  minor  principal  stresses,  respectively,  such  that  by 
selecting  the  values  of  k(3)  and  7(3)  become: 

k(3)  -  k  (3-25a) 

7(P)  -  tt/2  (3-25b) 


To  determine  the  sign  of  the  derivative  of  k  near  9  *  3>  Eqs.  (3-l6)  are  sub¬ 
stituted  into  Eqs.  (3-15)*  and  the  result  is  substituted  into  Eq.  (3-22): 


1-k  m  _ N  sin(3  -  29) _ 

1+k  cos  27(p  +  N(8in  3-29  cos  3)] 


(3-26) 


Taking  the  derivative  of  this  expression  with  respect  to  9,  noting 
Eq.  (3-25b),  and  that  N  according  to  Eq.  (3~l8)  does  not  depend  on  9,  one  finds 
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a  (i-k 

drp  u+k' 


4N2  sin  3  cos  g _ 

[p  +  N(sin  0-23  cos  3)]^ 


(3-27) 


cp«0 


N  being  a  real  number,  the  sign  of  the  above  expression  in  the  range  0  <  3  <  n 
is  necessarily  controlled  by  the  term  cos  3,  or 


sign  — 
dcp 


-  sign  [cos  3] 


(3-28) 


The  value  of  k(3)  *  k  is  by  definition  in  the  range  0  <  k  <  1,  as  consequence 

1-ic 

of  which  the  sign  of  the  derivative  of  k  and  of  -  are  opposite,  or  finally 

1+k 


sign  [k  (3)] 


sign  [-  cos  0] 


(3-29) 


Equation  (3-24)  requires  that  the  value  of  ,k(cp)  in  the  interior  of  the 
wedge  must  be  larger  than  the  value  at  the  boundary,  k(3)  *  k,  such  that  the 
derivative  k  (3)  must  be  negative.  Eq.  (3-29)  gives  therefore  cos  3  >  0, 
restricting  the  validity  of  the  elastic  solution  to  wedge  angles 

0  <  3  <  tt/2  (3-30) 


Next,  Eq.  (3-21)  defines  tan  (27)  as  function  of  the  stresses.  Using 
Eqs.  (3-15  and  l6)  one  finds 


tan  (2r)  . 


(3-3D 


In  the  range  of  3  given  by  Eq.  (3-30)  the  function  on  the  right-hand  side 
vanishes  for  cp  *  0  and  9  *  3>  assumes  negative  (positive)  values  for  0  <  cp  <  3/2 
(3/2  <  9  <  3)7  and  becomes  infinite  for  9  «  3/2.  Solving  for  27  gives  an  infi¬ 
nite  set  of  solutions.  Starting  with  a  value  27(0)  *  ±  nrr,  where  n  is  an  arbi¬ 
trary  integer,  2 7  decreases  from  above  starting  value  for  9  ■  0  to  an  end  value 
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27(3)  =  (±  n-1)  tt.  [It  is  easily  checked  that  tan  (27)  goes  through  the  seune 
range  of  values  as  the  right  side  of  Eq.  (3-31)0  Eq.  (3-25b),  defines  the 
value  7(3)  =  tt/2,  such  that  the  appropriate  value  of  n  must  be  n  =  2,  and 

7(0)  =  tt  (3-32) 


The  angle  7  defining  the  direction  of  the  principal  stress  decreases  gradually 
from  7(0)  =  tt  to  7(g)  -  tt/2. 

Equation  (3-32)  defines  the  position  of  the  principal  stress  for 

cp  =  0.  Referring  to  Fig.  15,  cr^  acts  in  the  radial  direction,  or  =  ar(0). 

The  other  principal  stress  is  then  =  0^(0).  These  relations,  together  with 
Eqs.  (3-15,  and  l8)  permit  determination  of  the  value  k(0)  from  Eq.  (3-23): 


k(0) 


CTa  yo) 


p  +  2N  sin  3 


(1-k)  +  2k 
(1-k)  +  (1+k) 


(3-33) 


The  angle  8  being  limited  by  Eq.  (3-30)  to  0  <  8  <  tt/2  and  the  material  pro¬ 
perty  k  being  limited,  0  <  k  <  1,  numerator  and  denominator  are  each  the  sum 
of  two  positive  terms.  Further,  in  the  permissible  range  of  8  the  value  of 
tan  8/8  >  0,  and  since  2k  <  X+k,  the  value  of  the  fraction  in  Eq.  (3-33)  cannot 
possibly  exceed  unity,  k(0)  <  1.  The  smallest  value  of  the  fraction  occurs 
when  tan  8/8  is  as  large  as  possible,  l.e.,  tan  8/8  -*  °°  ,  which  case  gives  a 
lower  limit,  k(0)  >  2k/(l+k).  In  turn,  1+k  <  2,  such  that  2k/(l+k)  >  k,  or 
finally 

1  >  k(0)  >  k  (3-3*0 

The  condition  for  elastic  behavior,  Eq.  ( 3-24 ),  is  therefore  satisfied  at  cp  =  0. 

As  next  step  the  sign  of  the  derivative  k ' ( 0 )  is  determined  by  differen¬ 
tiation  of  Eq.  (3-26).  Noting  Eq.  (3-32),  one  finds 
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—  (/il£N'  M  _  2p  N  cos  3 
dtp  1+k^  Q  [p  +  N  sin  9]* 


(3-35) 


In  the  range  (3-30),  the  values  of  cos  3  and  of  N,  the  latter  defined  by 

Eq.  (3-18),  are  positive,  such  that  the  above  derivative  is  negative.  The 

-  1-k 

sign  of  the  derivative  k  being  opposite  to  that  of  -  , 

i+ic 

sign  [k'(0)j  >  0  (3-36) 


The  results  obtained  so  fax  may  be  summarized  bv  considering  Fig.  1 6 
which  shows  a  typical  plot  of  k(cp).  The  end  points  of  the  curve  have  been 
determined  above,  k( 0 )  being  limited  by  the  inequality  (3-34),  while  k(9)  »  k. 
Equations  (3-29,  30  and  36)  define  the  tangents  to  the  curve  k(cp)  at  cp  »  0  and 
.  <p  *  3  in  the  manner  shown  in  Fig.  l6.  Between  these  end  points  k(cp)  is  defined 
'  by  Eq.  (3-26)  as  a  function  of  cp.  The  form  of  this  equation  is,  however,  not 
suitable  to  demonstrate  that  k(cp)  is  a  smooth  function,  as  indicated  in  Fig.  1 6. 
If  Eq.  (3-20)  is  used,  the  function  sin  (9  -  2tp)/  cos  2/  becomes  indefinite  for 
cp  *  9/2  when  7(9/2)  «  3n/4.  This  difficulty  is  overcome  by  computing  the  value 
,  of  sin  (9  -  2cp)  from  Eq.  (3-31)  and  substituting  the  result  into  Eq.  (3-26), 


1+k  sin  27  p  +  N(8in  9  -  2cp  cos  9) 

The  first  fraction  on  the  right  side  of  this  equation  is  obviously  smooth, 
because  7  varies  smoothly  from  7(0)  -  tt  to  7(0)  -  n/2,  see  Eq.  (3-32)  and  the 
adjoining  text.  Excluding  the  end  values  7(0)  and  7(9),  the  value  of  sin  27  <  0, 
and  because  N  >  0,  as  previously  noted  in  connection  with  Eq.  (3-35),  the  first 
fraction  is  necessarily  positive. 
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f 


f 

» 


To  recognize  the  character  of  the  second  fraction,  the  value  of  N, 
Eq.  (3-18),  is  substituted 


1 

p  +  N(sin  g  -  2cp  cos  3) 


23  cos  3  +  sin  3 
p[2(3-<p)  cos  3  +  sin  3J 


(3-38) 


In  the  range  0  <  3  <  tt/2  all  terms  on  the  right  side  are  positive,  and  the 
fraction  (3-38)  is  necessarily  positive.  Returning  to  Eq.  (3-37)  it  follows  that 


>  0  (3-39) 

1+k 


If  the  maximum  value  of  k,  indicated  in  Pig.  l6,  were  larger  than  unity,  the 

expression  would  become  negative,  such  that  Eq.  (3-39)  indicates  that 
1+k 


max  k(cp)  <  1 


(3-40) 


Allowing  for  the  minimum  of  ic(cp)  at  <p  ■  3,  min  k  -  k,  the  values  k(cp)  satisfy 
Eq.  (3-24).  It  has  therefore  been  proved  that  the  original  assumption  that 
slip  does  not  occur  in  the  wedges,  0  <  <p  <  3;  is  justified  and  furnishes  the 
solution  to  the  problem  considered. 


To  complete  the  analysis,  we  obtain  the  inclination  3  of  the  locking 

front  as  function  of  the  velocity  V  of  the  surface  pressure  from  Eqs.  (3-7  and  8) 

« 

used  in  the  elastic  case.  One  finds 


V2 


JL 

_ _ _ 

k  sin2  3  +  (l-k)  0  sin  3  cos  3 


(3-41) 
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Just  as  in  ti.e  elastic  case,  Eq.  (3“9)>  V  as  function  of  8  may  be  plotted, 

and  the  result,  shown  for  k  ■  l/3  in  Fig.  17,  is  quite  similar  to  the  one  for 

the  elastic  case,  except  that  Eq.  (3-30)  restricts  the  values  of  the  angle  3 

to  the  range  0  <  3  <  tt/2.  Real  values  of  V  are  found  everywhere  in  this  range. 

These  values  of  V  have  a  minimum  Vcr  ,  which  occurs  for  a  critical  wedge  angle, 

8  <  rr/2.  The  latter  may  be  found  from 

cr  ' 


tan  (28cr)  2(l-k) 

8  "  ‘  1+k 

cr 


(3-te) 


Figure  17  shows  that  for  values  of  V  near  Vcr  there  are  two  values  of  8, 
one  larger,  one  smaller  than  8  •  However,  using  the  reasoning  employed  in 

the  elastic  case  [in  the  text  following  Eq.  (3-9)] >  it  is  concluded  that  only 
wedge  angles  8  <  8  are  physically  meaningful.  The  inapplicable  portion  of 

vl 

the  curve  in  Fig.  17  is  shown  in  dashed  lines. 


The  fact  that  no  solution  has  been  obtained  if  V  <  V  is  again  ascribed 
to  the  fact  that  a  steady  state  solution  for  an  expanding  surface  load  need  not 
have  the  character  of  Fig.  11,  where  the  locking  front  and  the  face  of  the 
applied  pressure  move  together.  For  velocities  below  the  critical,  a  detached 
locking  front  forms,  Fig.  13,  and  the  steady  state  solution  is  given  by  [6] . 

In  applications  the  velocity  V  of  shock  waves  in  air  is  always  larger  than  V  , 
such  that  detached  locking  fronts  are  merely  an  academic  possibility.  The 
solutions  derived  above  apply  therefore  in  cases  of  practical  interest,  as  may 
be  seen  from  the  typical  example  at  the  end  of  Section  IV. 


It  is  necessary  to  remember  that  the  solution  just  obtained,  shown  in 
Fig.  17  for  k  «  l/3,  applies  only  if  the  inequality  (3-12)  is  satisfied.  For  the 
example,  k  *  l/3,  this  requires  v  <  l/4.  The  alternative  case,  when  the  inequa- 
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lity  (3-13)  applies  remains  to  be  discussed.  If  the  latter  inequality, 


applies,  the  stresses  at  and  immediately  behind  the  consolidation  front  «re 
elastic,  such  that  the  conditions  at  the  boundaries  of  the  wedge  are  exactly 
those  used  for  the  elastic  half -space,  Eqs.  (3-2)  and  (3-3)*  It  will  be 
shown  that  the  elastic  solution  found  in  Section  A  for  these  boundary  conditions, 
represents  a  state  of  stress  in  which  condition  (2-43)  is  satisfied  everywhere, 
i.e.  no  slip  occurs  in  any  location.  This  means,  that  the  solution  for  an  elas¬ 
tic  material  not  subject  to  slip  applies  in  case  of  the  inequality  (3-13)  also 
for  the  material  subject  to  slip. 

To  prove  above  contention,  it  is  noted  that  the  solution  in  Section  A 
(for  the  elastic  material)  may  be  obtained  from  the  solution  for  the  material 
subject  to  slip  by  replacing  the  value  k  by  the  expression  v/(l-v).  This  is 
due  to  the  fact  that  the  only  difference  in  the  formulation  of  the  two  solutions 
originates  from  the  boundary  conditions  for  ar(3)>  Eq.  (3-3a)  and  Eq.  (3-l4a), 
respectively,  where  in  one  case  the  term  v/(l-v)  occurs,  while  k  occurs  in  the 
other.  The  results  concerning  the  value  of  the  function  ic(cp)  obtained  above, 

Eqs.  (3-19)  to  (3-42),  can  therefore  be  utilized  without  further  computations. 
Figure  18  shows  a  typical  distribution  of  k(cp).  At  the  boundary  cp  «  3,  the 
ratio  k(3)  of  the  two  principal  stresses  is,  by  virtue  of  Eq.  (3~3a), 

l 

k(3)-i^;>k  (3-43) 

The  end  point  of  the  k-curve  at  3  is  therefore  above  k.  Further,  the  value 
k(3)  must  be  less  than  unity,  because  Poisson's  ratio  is  always  less  than  l/2. 


Replacing  k  by  v/(l-v)  does  not  affect  the  results  obtained  in  Eq.  (3-29 )> 
Eq .  (3-3*0*  Eq*  (3-36)  and  Eq.  (3-40),  defining,  respectively,  ic  (8),  k(0),  k  (0) 
and  max  k.  The  values  k (cp )  lie  therefore  everywhere  between  the  limits 

k  <  k(cp)  <  1  (3-44) 

confirming  the  claim  made  above  that  slip  does  not  occur  anywhere. 


While  this  was  previously  stated,  it  is  stressed  again  that  all  results 
in  this  Section,  for  Case  A  as  well  as  Case  B,  apply  only  for  velocities  V  below 
a  certain  limit.  This  is  due  to  the  fact  that,  in  similarity  to  [1],  the 
equations  of  motions  in  the  elastic,  wedge-shaped  regions  were  approximated  by 
Eqs.  (3-1)*  This  approximation  is  reasonable  only  if  the  velocity  V  is  a 
fraction  of  the  velocity  of  P-waves 


2  E(l-v) 

Cp  *  p(l+v)(l-2v) 


(3-45) 


where  E  and  v  are  the  elastic  constants  describing  the  behavior  of  the  material 

after  compaction,  if  no  slip  occurs.  It  was  shown  in  [1]  that  the  approximate 

Eqs.  (3-1)  may  be  used  if  V/cp  <  0.20,  such  that  the  present  results  will  apply 

in  this  range.  However,  for  values  of  the  parameters  of  interest  in  applications, 

the  locking  front  makes  a  very  small  angle  with  the  surface,  0  <  0.20  (see  the 

example  at  the  end  of  Section  IV).  A  comparison  with  a  refined  solution  shows 

that  for  such  small  angles  of  8,  and  for  reasonable  values  of  Poisson’s  ratio, 

v  <  1/3,  the  use  of  Eqs.  (3-1)  is  permissible  even  for  values  of  v/c  higher 

P 

than  0.20.  The  present  solution  is  expected  to  give  good  results  if  the  in¬ 
equalities 

v  <  1/3  ,  8  <  0.20  (3-46) 

and 

V/cp  <  0.50  (3-47) 
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are  satisfied.  These  inequalities  insure  that  V  >  cg  ,  where  cg  ■  /G /p  is  the 

velocity  of  shear  waves  in  the  elastic  (compacted)  material.  cg  being  about 

0.60  to  0.70  c  ,  the  present  solution  applies  for  a  major  portion  of  the 
P 

"subseisraic"  range,  V  <  cg  . 
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IV  DISCUSSION  OF  RESULTS 


In  the  previous  Section,  Case  B,  the  stresses  in  a  half-space  were 
obtained  for  locking  materials  of  elastic -Coulomb  slip  behavior  after  com¬ 
paction.  It  was  found  that  two  situations  must  be  distinguished,  and  the 
results  differ,  depending  on  whether 


1-v 


< 


k 


or 


v 

1-v 


> 


k 


(4-1) 


(4-2) 


The  solutions  obtained  indicate  that  a  locking  front  occurs,  inclined 

at  an  angle  0,  Fig,  19,  when  the  velocity  V  of  the  applied  surface  load  p 

exceeds  a  critical  value  V  .  The  value  o  (0)  of  the  normal  stress  at  the 

ci*  cp 

front,  the  relations  between  V  and  0,  and  the  value  of  0  are  given,  respectively, 
by  Eqs,  (3-17;  4l,  42)  in  the  case  of  the  inequality  (4-1),  and  Eqs„  (3-6,  9,  11) 
in  the  case  of  (4-2).  It  was  further  concluded  that  only  values  of  0  <  0cr  are 
physically  meaningful.  Figure  17  shows  a  typical  plot  of  V(0)  for  k  =  l/3 
when  Eq.  (4-1)  applies,  while  Figs.  10  and  l4  present  V(0)  and  0^(0),  respectively, 
for  v  =  l/3  in  case  Eq.  (4-2)  holds. 


The  details  of  the  solution  show  that,  regardless  whether  Eq.  (4-1)  or 
(4-2)  applies,  no  slip  occurs  in  the  interior  of  the  wedge-shaped  areas 
0  <  cp  <  0,  which  act  entirely  elastic.  However,  in  the  case  of  Eq.  (4-1)  slip  does 
occur  within  the  (theoretically  infinitely  thin)  compaction  front.  In  the  case  of 
Eq.  (4-2)  slip  does  not  occur  at  all,  not  even  at  this  front. 
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From  the  results  concerning  the  stress  field,  other  quantities  of 
interest  can  be  obtained.  One  important  quantity  is  the  permanent  strain 
produced  at  any  point  due  to  the  passing  of  the  pressure  front  on  the 
surface.  The  significant  portion  of  the  permanent  strain  occurs  during  the 
passing  of  the  locking  front,  the  principal  strains  being 

«1  =  “  «  i  «2  =  °  (4-3) 

where  c  is  the  locking  strain,  while  are  inclined  at  angles  to  the 

surface  of  90°  -  3,  and  3>  respectively.  Changes  of  strain  subsequent  to 
tne  passing  of  the  locking  front  are  negligible  compared  to  the  locking  strain  e, 
because  in  an  ideal  locking  material,  when  no  slip  occurs,  the  post-locking 
elastic  strains  are  of  the  order  |p/E|  <  <  c.  Equations  (4-3)  define  there¬ 
fore  the  total  permanent  strain.  Without  interaction  of  structure  and  free 
field  an  originally  circular  structure,  Fig.  19,  would  be  deformed  to  an 
ellipse  of  the  axis  ratio  1  -  e,  and  the  original  volume  of  the  structure 
would  be  reduced  by  a  factor  1  -  *.  A  cylindrical  shell,  e.g.,  will  not 
strongly  resist  the  deformation  into  an  ellipse,  but  may  be  able  to  resist  a 
volume  change  materially.  Predictions  on  the  behavior  of  structures  in  locking 
materials  require  therefore  further  analysis  of  interaction  effects  which  will 
modify  the  free  field  results  on  stress  and  strain  obtained  above. 

One  may  also  obtain  the  accelerations  a  and  particle  velocities  u  which 
correspond  to  the  stress  field  obtained  in  Section  III.  Again,  because  we  con¬ 
sider  an  ideal  locking  material,  the  significant  effect  is  the  change  in  particle 
velocity  at  the  locking  front.  Ahead  of  the  front,  u  *  0.  To  obtain  the  particle 
velocity  u^  behind  the  locking  front,  it  is  noted  that  the  front  is  plane,  such 
that  the  well  known  relation  applying  for  one  dimensional  wave  propagation 
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UL  *  eU 


(4-1.) 


may  be  used,  where  (J  is  the  velocity  of  the  front,  and  uT  is  the  particle 

li 

velocity  at  right  angles  to  the  front  after  locking.  The  velocity  U  is 
given  by  Eq.  (2-39 )>  where  a  means  the  normal  stress  at  the  front,  in  the 

aA 

present  situation  o  s  a  (0).  The  particle  velocity  uT  is  therefore 

XX  CD  J_i 


-  «oJ3) 


(4-5) 


where  0^(3)  is  given  by  the  appropriate  Eqs.  (3-6)  or  (3-17).  Fig.  20  gives 
the  time  history  of  the  particle  velocity  u(t),  its  direction  making  the  angle 
90°  -  3  with  the  surface.  Thf  jump  occurs  when  the  locking  front  passes  the 
point  of  observation.  The  accelerat_on  a  vanishes  except  at  the  instance  of 
passing  of  the  front,  where  a  *  «  ,  hut  the  order  of  the  infinity  is  such  that 


\a  dt  *  uT 


(4-6) 


The  direction  of  a  is  the  same  as  that  of  u(t).  This  result  may  be  used  for 
the  determination  of  shock  factors. 


For  materials  which  behave  elastically  after  locking,  no  separate  dis¬ 
cussion  is  necessary.  Statements  and  results  apply,  as  given  for  materials 
with  '.lip  when  the  inequality  (4-2)  holds. 


Consider  a  soil  of  density  p  *  1.8  x  10~^  lb.sec?/ln^  at  a  pressure  level 

p 

p  a  200  lb. /in.,  when  V  »  3970  ft. /sec.  Assume  a  material  with  slip,  k  «  1/3, 
v  *  l/5,  e  *  0.02.  The  inequality  (4-1)  holds,  such  that  Eqs.  (3-17;  4l,  42) 
apply. 
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To  ascertain  that  V  <  V,  one  finds  from  Eq.  (3-42):  B  *  1.1^5  rad., 

cr  cr 

and  from  Eq.  (3-4l)  the  corresponding  value  Vcr  »  820  ft. /sec.  <  V.  Eq.  (3-41) 

will  therefore  have  a  root  B  <  B  for  the  above  parameters  and  the  analysis 

cr 

presented  above  furnishes  the  stresses,  velocities,  etc.,  in  this  situation. 
Solving  Eq.  (3-41 )  gives  B  *  0.070  rad.,  and  Eq.  (3-17)  furnishes  the 

o 

pressure  at  the  shock  front,  a  (B)  ■  201  lb. /in.  This  pressure  differs  only 
very  slightly  from  the  surface  pressure  p.  This  is  due  to  the  fact  that  B  is 
quite  small  whenever  V  is  very  much  larger  than  Vcr  .  The  particle  velocity 
behind  the  locking  front,  found  from  Eq.  (4-5)>  is  uL  «  20  ft. /sec. 

Due  to  the  fact  that  the  angle  B  is  small,  B  <  <  1>  the  solution  differs 

only  little  from  the  one  dimensional  one.  In  the  latter  the  stress  at  the 

2 

locking  front  for  an  applied  step  pressure  p  *  200  lb. /in.  would  also  be 
200  ib./inT,  while  the  two  dimensional  solution  gives  or  (P)  *  201  lb. /in. 

The  analysis  presented  in  Section  III  is  restricted  to  values  of  the 
parameters  which  satisfy  the  inequalities  (3-46,  47 )•  Eqs.  (3-46)  ob” 
viousiy  satisfied.  Introducing  Eq.  (3-45)  into  (3-47)  leads  to  the  requirement 

E  >  v2  (4-7) 

6  2 

which  gives  in  the  present  case  E  >  1.5  x  10w  lb. /in.  The  value  E  in  this 
inequality  is  Young's  modulus  in  the  contacted  material.  At  the  present  time 
there  is  little  factual  information  on  this  value  of  E,  for  loading  or  'unloading, 
subsequent  to  compaction. 


In  the  example  discussed  above,  the  velocity  was  substantially 
smaller  than  the  shock  velocity  V  in  the  air.  This  situation  holds  generally 
for  a  wide  rang*,  of  t  and  p.  This  is  demonstrated  by  the  following  table: 


Table  of  Vcr  as  function  of  p  and  e  for  k  =  1/3.  v  <  l/4 


c 

mm 

p  =  500  lb ./in? 

0.01 

820  ft. /sec. 

1840  ft. /sec. 

0.10 

260  ft. /sec. 

580  ft. /sec. 

The  values  of  V 

cr 

and  6130  ft. /sec. 


are  much  smaller  than  the  shock  velocities  in  air,  V  ~  2920 

p 

for  p  =  100  and  500  lb. /in.,  respectively. 


V  SUMMARY  OF  CONCLUSIONS 


As  an  introductory  step,  Section  II  contains  a  study  of  the  conditions 
on  locking  fronts  in  problems  of  two  and  three  dimensional  wave  propagation. 
Considering  a  locking  material  as  a  limiting  case  of  a  nonlinear  (hardening) 
elastic  material,  it  was  found  that  the  normal  stress  at  a  locking  front  must 
be  a  principal  stress,  such  that  shear  stresses  vanish. 


Based  on  the  study  of  the  conditions  on  locking  fronts,  the  paper  con¬ 
tains  an  analysis,  giving  stresses,  strains  and  particle  velocities  due  to  a 
progressing  step  wave  on  the  surface  of  a  half-space  of  a  locking  material, 
which  after  compaction  is  elastic  but  subject  to  slip.  The  slip  conditions, 
previously  derived  in  ( 1 J ,  define  the  states  of  stress  necessary  to  overcome 
internal  Coulomb  friction  in  terms  of  a  parameter  k,  which  is  the  ratio  of  the 
minor  (compressive)  principal  stress  divided  by  the  major  one. 


The  details  of  the  results,  given  in  Section  IV,  depend  on  the  relative 
values  of  k,  and  of  Poisson's  ratio  v  as  expressed  by  the  inequalities  (4-1) 
and  (4-2).  There  is  a  locking  front,  inclined  at  an  angle  0,  Fig.  19, 
followed  by  a  gradually  varying  stress  field.  Slip  occurs  only  at  the  locking 
front  if  Eq.  (4-1)  applies,  and  not  at  all  if  Eq.  (4-2)  applies. 


Because  of  a  simplification  in  the  analysis,  the  results  are  restricted 

with  respect  to  the  velocity  V  of  the  surface  pressure.  In  the  case  of  air- 

blast  on  the  surface  of  a  granular  soil,  the  "Mach  number"  must  be  limited 

to  V/^  <  0.50,  where  c^  is  the  velocity  of  P-waves  in  the  compacted  material. 

The  velocity  c  of  shea"  waves  being  only  slightly  larger  than  0.5  c  ,  the 
s  p 


range  of  validity  of  the  present  solution  covers  most  of  the  subseismic  range, 

V  <  c  .  Solutions  for  higher  values  of  V,  particularly  for  the  superseismic 
s 

range,  V  >  c  ,  require  a  different  approach,  left  for  a  future  investigation. 

P 

As  a  byproduct  of  the  analysis,  the  solution  for  a  locking  material 
which,  after  compaction,  acts  elastically,  without  slip,  has  also  been  obtained, 
subject  to  the  limitation  V/c^  <  0*50. 

The  example  at  the  end  of  Section  IV  indicates  that,  for  airblast,  the 
results  of  the  two  dimensional  analysis  differ  only  minutely  from  the  results 
of  the  much  simpler  one  dimensional  analysis.  One  can  expect  that  this  con¬ 
clusion  will  also  apply  for  a  decaying  surface  pressure,  and  one  dimensional 
results  may  therefore  be  used  approximately  for  two  dimensional  problems. 


k3 


FRONT  AT  TIME  t 


FRONT  AT  TIME  t  +  dt 


-cr 
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